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Abstract 

Imaging of sources from data within the intensity interferometry is discussed. 
In the two-pion case, the relative pion source function may be determined by 
Fourier transforming the correlation function. In the proton-proton case, 

the discretized source function may be fitted to the correlation data. 
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Phase interferometry [I]] is capable of delivering star images in astronomy, e.g. that of 
Betelgeuse 0. Intensity interferometry, as is applied in nuclear physics, has primarily been 
used to determine radii of stars in astronomy and of particle emitting regions in heavy-ion 
reactions. In heavy-ion reactions, this determination has most often been done by fitting the 


low-momentum two-particle correlation functions under the assumption of Gaussian-shaped 
emitting regions . Similarly, source lifetimes have been inferred by considering Gaussian- 


shaped distributions of emission times. However, it was not always clear that the analyzed 


data truly narrows down the lifet im e for short-lived sources. Bertsch [d| and Mrowczynski || 


noted that an integral over the low-momentum correlation function can yield the value of the 
source function at the zero relative distance. Beyond this, no attempts were made to image 
the source function from the reaction data. Typically, comparison to reaction simulations is 
carried out by generating the correlation functions from simulated events. 

In this letter, we investigate the feasibility of direct imaging of the source from reaction 
data, within the intensity interferometry. Given a two-particle correlation function, this 
represents an inversion problem. We consider two examples: that of like-charged-pion and 
that of proton-proton interferometry. In the like-pion case, the relative source function 
is a Fourier transformation of the two-particle correlation function. In the proton case, 
the procedure is more involved. In this letter, we first discuss the relation of the correlation 
function to the source function, then the imaging, and finally some information contained 
in the images. 

Under the assumption of the approximate independence of processes leading to two 
particles in the final state and the assumption of the weak dependence of the product of 
single-particle sources D on momenta important for correlations, the two-particle correlation 
function may be represented as [Sj,j7| 



( 1 ) 


Here Sp (r) is the distribution of relative separation of emission points for the two particles, 
in their center of mass. In terms of single-particle sources, 
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<Sp(r) = / dUdti dt 2 D{ 0, R + r/2, D{ 0, R — r/2, t 2 ) 


( 2 ) 


The momentum P above is the total momentum of the pair and q is the cm relative momen¬ 
tum. The integrations in (|1) and in (2) are over cm variables; the single-particles sources, 
D = D / J dr dt D , are taken in the pair cm frame. The function D is the distribution of 
last collisions for an emitted particle in space, time, and momentum. Both S and D are 
normalized to 1. For Klein-Gordon fields, with (□ + m^)(j){x) = —j(x), D may be written 
in terms of single-particle self-energies as 

1 


D{ p, r, t) = —- n<(p, E p , r, t) exp 

Ztip 


2 E p Jt 


POO 

J dt' (—2)Im fl + (p, E p , r + v p (f' - t),t') 


( 3 ) 


where ^{x.x') = (j(x')j(x)) irre d, and (-2)ImfI + (a;, x') = ([j(x),j(x , )]) i „ ed . For the 
Schrodinger fields, with (i-j^ + T(x) = j(x), the analogous result is 


POO 

D(p, r, t) = =FfE < (p, E p , r,t) exp - J dt'T (p, E p , r + v p (f' — t),t') 


( 4 ) 


where =!=*£< is the single-particle production rate, =Ff£ < (a:, x') = (j(x')j(x)) i rre d , and V is 
the damping rate. For particles with spin, the modulus of the wavefunction in ([[]) should 
be averaged over spin directions. For low r within the source, deviations from (^) could be 
expected in the presence of short-range repulsion. Results (JT|) , (|3|), and (|4|) ignore the final- 
state refraction. We discuss the effects of the Coulomb field of a source in a separate paper || . 
For pion-pion correlations, both these effects and the effects of the pion-pion Coulomb 
interaction may be approximately removed directly from the correlation function ||. 

While the correlation function C has an indirect dependence on other quantities, it is 
shaped directly by S in the reaction, i.e. by the relative distribution of emission points 
for two particles with similar momenta, in their center of mass. For like particles, S is a 
symmetric function. This is not generally the case for distinct particles. Depending on the 
circumstances in a reaction, S may range from isotropic, for prompt emission, to strongly 
elongated along P, and possibly even bone-shaped, for emission from a long-lived source. In 
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this last case, the transverse size would be close to the size of a nucleus at short distances 
along P and widen at larger distances due to the zigzaging of the source under the recoil 
caused by emission. The lifet im e might be read off from the elongation of S. 

With (I|), the goal of the imaging is the determination of S given C. Given that the inter¬ 
esting part of C is its deviation from 1, we may subtract 1 from both sides of (lj) obtaining 

Cp(q) - 1 = J dr (|3>q~ ) (r)| 2 - l) S'p(r) = J dr K(q, r) 5 P (r), (5) 

where K = | < h^’| 2 — 1. The problem of imaging then reduces to the inversion of K. A dif¬ 
ficulty may arise from the presence of a kernel (or null-space) of K, i.e. the subspace of 
functions that the operator K turns to zero. The projection of S onto the kernel cannot 
be restored within imaging. It will become apparent that, in the case of like particles, the 
kernel of K is empty. For unlike particles, when one of the particles is neutral, the imaging 
may not be able to restore portions of S for large particle separations as K approaches zero 
for large separations; a particular severe situation occurs for 1$^'| 2 — 1, when the whole 
space of functions becomes the kernel. 

For like-pion pairs, Eq. ([5]) may be written as 

Cp(q) — 1 = / dr cos (2q ■ r) Sp(r), (6) 

if we ignore interactions between the pions. Given that S is symmetric, the Fourier transform 
can be inverted to yield 

-Sp(r) = j dq cos (2q ■ r) (C P (q) - 1) . (7) 

The directions that we use in the analysis of S are, in the system frame, outward along the 
transverse momentum of the pair, longitudinal along the beam, and the remaining direction, 
termed transverse. With equation (7]), we can find the angular moments of the source and the 
correlation function. If we introduce C(q) = C Xm (q) Y Am (f] q ) and an analogous 

representation for S, then we get the relation 

s£ m (r) = IzAdi r dqq 2 Ja (2 qr) (C£ m (r) - S ml> ) . (8) 

7 r- Jo 
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Due to the symmetry of S and C, only even A appear in the angular expansion of these 
functions. Since both functions are real, the moments satisfy (C Am )* = (—l) m C x ~ m . The re¬ 
lation (^]) between the angular moments may help in analyzing the three-dimensional data. 
In particular, this relation shows that the angle-averaged correlation function C 00 (q) = C(q) 
reflects the angle-averaged source S 00 (r ) = S(r), 


2 r°° 

r Sp(r) = — / dq q sin (2qr) (Cp(q) — 1). 
tt z Jo 


(9) 


As a specific example of the source extraction, in Fig. [I] we present the relative angle- 
averaged 7source-function. It is determined by applying Eq. H) to the data of Ref. El 


for central 10.8 GeV/c Au + Au. Prior to the Fourier transformation in (J|), the data were 
corrected for the Coulomb interaction between the two pions and between the pions and 
the source 0. The integration in (||) for Fig. |I| was cut off at g max ~ 50 MeV/c, giving 
a resolution in the relative distance in the figure of Ar rv./ 1/2 g max rv./ 2.0 fm. The largest 
r that may be considered is 1/2 A q, where A q is the momentum resolution of the data 
(A q = 5 MeV/c in the case of |10|). 

In the general case, the respective angular moments of C and S are also directly related. 
The spin-averaged operator, K, only depends on the angle between q and r, and not on 
their separate directions. Thus, the averaged K may be expanded: K( q, r) = Xfy(2A + 
1) K\(q,r ) P x (cos6). A relation between the moments follows from (0), 


roo 

C Xm {q ) - <5 ao 5 m0 = 4tt / dr r 2 K x (q , r) S Xm (r) 

Jo 


( 10 ) 


In the like-nucleon and like-charged-pion cases, the A = 0 operator is 


K o(q, = ( 2 J + !) [pf a ( r )) 2 - 1 1 


( 11 ) 


jsW 


where g^ s is the radial wave function with outgoing asymptotic angular momentum l. If the 
correlations are of purely Coulomb origin, such as between intermediate-mass fragments ]TI 


the operator is K 0 (q,r) = 6{r — r c ) (1 — rjr) 1 ^ 2 — 1 in the classical limit. Here r c is the 
distance of closest approach, r c (q) = 2/iZ\ Z 2 e 2 /g 2 . 
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We determine the source by discretizing the functions and integrals in equations such as 
(0). ©, or ©> and fitting the discretized values of S. We illustrate this by analyzing 
the proton-proton correlation data 0,|F|] from the 75 MeV/nucleon 14 N + 2 'A1 reaction, 
displayed in Fig. ^]. 

With the data [|T3J] averaged over the q directions, we concentrate on the relation © 
between the angle-averaged S and C. At high relative-momenta, the assumptions leading 
to ( 0 ) may break down. Further, at high momentum correlations within the source, such as 
associated with the ansisotropic flow [|Fj, may play a role. In addition, the event selection 


for singles in [13] becomes an issue. Thus we restrict the region for source determination 


to less than g max ~ 80 MeV/c. Given the issues of experimental resolution |T3|, we restrict 
the region of analysis from below to q > q mir , = 10 MeV/c. This sets a lower limit on the 
resolution within the source of 1/2(g max — Q m\ n ) — 1.4 fin. Therefore, we settle with a deter¬ 
mination of the source values at points separated by a coarser A r = 1.8 fin. We represent 
the source function as 


Sp(r) - s p( r k) g(r - r k ), 


( 12 ) 


where r k — (k — |)A r and g is a profile function, g{x) = 1 for |x| < Ar/2, and g{x) = 0 
otherwise. Then from ([lOi), we get 


Cp(q) - 1 ~ S P (r k ), 


(13) 


where w k = 47 rf drr 2 K 0 (q,r) g(r — r k ). We determine the pp wavefunctions for £, i' < 2 
from the Schrodinger equation with the regularized Reid soft-core potential |R|. Finally we 
minimize x‘p = Sj((C , p XP (Q , j) ~ Cp(qj))/(Jj ) 2 by varying S(r k ), subject to the conditions that 
S(r k ) > 0 and that S is normalized to 1. Initially, we aimed to determine the source up to 
t max ~ 30 fin, but we found that the fits favored S consistent with zero at higher r. Thus, 
we were able to reduce r max to 16.2 fm without an appreciable worsening in the fits. 

The source functions extracted from the data 00 are shown for the three total- 
momentum gates in Fig. |3] together with the source functions determined directly within 
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Boltzmann-equation reaction-simulations ||T7>| . p~2| . Unlike what was done in [12[], we initialize 
the nuclear matter density in the model by solving the Thomas-Fermi equations. The errors 
on the extracted S include the uncertainty from varying g max in the vicinity of 80 MeV/c. 
The values of y 2 per degree of freedom are 1.1 for the highest momentum gate, 1.8 for the 
lowest, and 4.9 for the intermediate one with the lowest errors on the data. 

We now discuss the relative proton source functions extracted from the data and calcu¬ 
lated in the model, together with the information contained in source functions. The change 
in the relative distribution of emission points in Fig. from a compact form at high pro¬ 
ton momenta to an extended form at low momenta, demonstrate the presence of space- 
momentum correlations within the 14 N + 2 'A1 reaction. At intermediate and low momenta, 
neither the relative proton distributions from the data, nor those from the model, can be 
well approximated by Gaussians. This is in contrast to the distribution in Fig. [lj. and to 
a degree, to distributions for high momenta in Fig. [3]. Overall, the transport model yields 
distributions in a near-quantitative agreement with the experimental sources. The model 
source seems to systematically underestimate the values from the data only at the shortest r. 

For a rapid freeze-out, the single-particle source is given by D{ p, r, t ) ~ /(p, r) 5(t — t 0 ) 
where / is the Wigner function. Assuming weak directional correlations between the pair to¬ 
tal and relative momenta, and between the spatial and momentum variables, the momentum 
average of S approximates the relative distribution of emission points for any two particles 
from the reaction, and not just for the particles with close momenta. Thus given a rapid 
freeze-out, the relative distribution for any two particles is 

f dP dp dR /(P/2 + p, R + r/2) /(P/2 - p, R - r/2) 


5(r) = 

/ dpi dri/(pi, ri) f dp 2 dr 2 /(p 2 , r 2 ) 

We now rewrite and expand the expression in the numerator in (|l4]) : 

/(P/2 + p, R + r/2) /(P/2 - p, R - r/2) 

N /(P/2 + P, R + r/2) /(P/2 - p, R - r/2) 


(14) 


= / dri/(pi,ri) / dr 2 f(p 2 ,r 2 


= j dr[ /(p 1} ri) J dr 2 /(p 2 , r' 2 ) I 1 + p — + ... 


f dr[ /(P/2 + p, ri) / dr 2 /(P/2 - p, r') 
d 
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K /(P/2 + p', R + r/2) /(P/2 - p ; , R - r/2) 

/ dri /(P/2 + p', ri) / dr' 2 /(P/2 - p', r' 2 ) 

The gradient term must be proportional to a combination of the vectors P, rj, and r 2 . 
For weak directional correlations, it then averages to zero under the integration in (ED- 
Inserting (ED into (0 and keeping the leading term, we obtain 


( 15 ) 


, 1 f , dN dN . . .. 

o(r) - TT5. I “Pi“P2 -JZT XT 7P*5p(r + n P ( 7 p - l)(n P r)), 


(16) 


N 2 J ,ri c/pi dp 2 

where N is particle multiplicity and n P = P/P. The argument of Sp has been written in 
the cm frame of an emitted pair and y P is the Lorentz factor for the transformation from 
the system frame to the pair cm frame. 

Generally the relative distribution of emission points for any two particles at r —> 0 
gives an average freeze-out density when multiplied by N — 1. If the assumptions above are 
valid, this density may be obtained by multiplying the average (|I6|) of Sp(r —> 0) by N — 1. 


The transport calculations fllfj] indicate that the measured |F| coincidence cross sections 
for the 14 N + 2 'A1 reaction are dominated by rather central collisions with b ~ 2.8 fm. 
The chance of detecting two particles at wide angles simultaneously is large only for these 
collisions. The rms nucleon cm momentum in central collisions is ~ 185 MeV/c and at 25° 
this corresponds to a lab momentum of ~ 320 MeV/c for the nucleon or ~ 640 MeV/c 
total for the pair. Thus, the results for the intermediate-momentum gate in Fig. |3] best 
represent the average situation in the central collision. Presuming that the relative spatial 
distributions of other particles to protons is similiar to that of two protons, we arrive at 
an average nuclear density in the vicinity of any emitted proton of 17 x 5(r —> 0) ~ 
17 x 0.0015 fm~ 3 = 0.025 fm -3 = 0.16 no- Here we assume the participants to have a 
total mass of 18 following from the fireball geometry at 6 ~ 2.8 fm. The directional space- 
momentum correlations anticipated in collisions due to a collective motion, to shadowing, 
or to emission that is most likely not instantaneous, make this value an upper limit on the 
freeze-out density. 

Irrespective of any correlations or of the validity of instantaneous freeze-out, the product 
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of the r —> 0 source-function and the momentum distribution yields the space average of the 
phase-space occupancy at freeze-out, 


, t ( ^\ ( 2?r ) 3 E p dN c ( n\ 

</(p)>= 2JTT m^ S2 » (r ^ 0) - 


(17) 


(see also ||). Refs. 00 only give the inclusive proton cross-sections in the 14 N + 2 ' A1 
reaction at two angles. Furthermore, the cross sections include large contributions from 
peripheral events. Under these circumstances, we instead use the thermal distribution 


diV t h/dp oc 1 j{z 


-1 pP 2 /2mT 


+ 1) for the central events in formula (p~7[). Here, z is set from 


the requirement of maximum entropy. For ~ 9 participant protons at b = 2.8 fm in the 
14 N + 2 ' A1 reaction, this requirement gives z ~ 1.10 and T ps 10.2 MeV. Equation (|T7|) 
can now be used to determine the phase-space average of the occupancy at freeze-out, 
(/) — / dp ((/(p))) 2 / / dp (/(p)), and to estimate the entropy per nucleon, 


S_ _ I dp ((/( p)) log ((/( p))) - (1 - </(p))) log (1 - (/(P)))) 
A ~ /dp(/(p)> 


Use of the thermal momentum distribution yields (/) ~ 0.23 and S/A ~ 2.7 for the 14 N + 
2 ' A1 reaction. For a distribution with nonequilibrium features, these values should represent 
the lower limit on the average occupation and the upper limit on the entropy. Indeed, for 
the source from the transport model and a thermal distribution, we find an entropy about 
0.5 per nucleon higher than the entropy calculated directly within the model. 

We have shown how to determine the relative source functions for particles directly from 
correlation data. When going beyond the Gaussian fitting, direct Fourier inversion, or simply 
fitting the source, the singular value decomposition ]17| might be useful. This method allows 
one to determine what kind of source parametrization can be narrowed practically down by 
the correlation data. 
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FIGURES 

FIG. 1. Relative distribution of emission points for negative pion determined from the corre¬ 
lation data of Ref. juj] . Prior to the integration in Eq. © , the relative momenta in the correlation 
function have been increased by 8% to correct for the Coulomb effect of the source; similar reduction 
of the momenta in the positive-pion correlation function brings the two correlation functions (or the 
source functions) into rough agreement with each other. The shaded area represents uncertainty 
in the source function associated with the uncertainty in the determination of the ir~ correlation 
function © and due to the choice of q ma v ~ 50 MeV/c. 


FIG. 2. Two-proton correlation function for the 14 N + 2 ' A1 reaction at 75 MeV/nucleon. 


The symbols represent data [ 13,12) for three gates of total momentum imposed on protons emitted 
in the vicinity of 6 \ a b = 25°. The lines represent the correlation function for the extracted sources 
displayed in Fig. 0. 


FIG. 3. Relative source function for protons emitted from the 14 N 
in the vicinity of @i a b = 25°, within three total momentum intervals, 
function extracted from the data Jk|. 12]. Open circles represent the 
the Boltzmann-equation model |1(|. 


+ 2 'A1 at 75 MeV/nucleon, 
Filled circles represent the 
function determined within 
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